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Higher Loop Results for the Plaquette, Using the Clover and Overlap 
Actions* 

A. Athenodorou a , H. Panagopoulos a , A. Tsapalis a 

a Department of Physics, University of Cyprus, Nicosia CY-1678, Cyprus 

We calculate the perturbative value of the free energy in QCD on the lattice. This quantity is directly related 
to the average plaquette. 

Our calculation is done to 3 loops using the clover action for fermions; the results are presented for arbitrary 
values of the clover coefficient, and for a wide range of fermion masses. 

In addition, we calculate the 2 loop result for the same quantity, using the overlap action. 



1. INTRODUCTION 

In this work, we compute the perturbative ex- 
pansion of the average plaquette, in SU(N) gauge 
theory with Nf fermion flavours. We present 
separate calculations using the clover action (3- 
loops), and the overlap action (2-loops) pQ. The 
simpler case of Wilson fermions was performed in 

m 

The average plaquette can be related to the 
perturbative free energy of lattice QCD, as well 
as to the expectation value of the action. 

The results can be used: a) In improved scal- 
ing schemes, using an appropriately defined effec- 
tive coupling, b) In long standing efforts, starting 
with P] , to determine the value of the gluon con- 
densate, c) In studies of the intcrquark potential 
0| . d) As a test of perturbation theory, at its 
limits of applicability. 

In standard notation, the action consists of a 
gluonic and a fermionic part:: 



S — Sq +Sf, 

s G = pY. E ^ p ) 

p 

with 

E G (P) = 1-Re Tr(P)/A 



(1) 



(2) 



The sum runs over all lxl plaquettes P. The 
fermionic action Sf, in the clover case, contains 
the Wilson term with bare fermionic mass m, and 



the standard clover term multiplied by the coeffi- 
cient csw) which is a free parameter in the present 
work; it is normally tuned in a way as to minimize 
O(a) effects. 

The average value of the action density, S/V, is 
directly related to the average plaquette; in par- 
ticular, for the gluonic part we have: 



(S G /V) = 6f3(E G (P)). 



(3) 



As for (Sf/V), it is trivial in any action which is 
bilinear in the fermion fields 0, and leads to: 

(S F /V) = -ANN f (4) 

We will calculate (E G ) in perturbation theory: 

(E G ) = Cl g 2 + c 2 g 4 + c 3 g e + --- (5) 

The n-loop coefficient can be written as c„ = 
c n + c n where c„ is the contribution of diagrams 
without fermion loops and comes from dia- 
grams containing fermions. The coefficients c„ 
have been known for some time up to 3 loops |5I2| . 
The coefficients c„ are also known to 3 loops for 
Wilson fermions [2]; in the present work we ex- 
tend this computation to clover fermions. 

The calculation of c n proceeds most con- 
veniently by computing first the free energy 
— (lnZ)/V, where Z is the full partition function 
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Z = / [DUViPiViPi] exp(-S). 



(Eg) 



(6) 



The average of E G is then extracted as follows 
1 d flnZ 



(7) 
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In particular, the perturbative expansion of 
{\nZ)/V : 

(ln^ = do-^^ln/3 + | + | + -(8) 

leads immediately to the relations: 
c 2 = d 1 /{2AN 2 ) 1 c 3 = d 2 /(24N 3 ). 

2. CLOVER FERMIONS 

Up to 3 loops, a total of 62 diagrams con- 
tribute, of which 26 involve fermionic loops. The 
involved algebra of the lattice perturbation the- 
ory was carried out using our computer package 
in Mathematica. The value for each diagram is 
computed numerically for a sequence of finite lat- 
tice sizes. Diagrams must be grouped in several 
infrared-finite sets, before extrapolating their val- 
ues to infinite lattice size; extrapolation leads to a 
(small) systematic error, which is estimated quite 
accurately 

The pure gluonic contributions are already 
known: 



N 2 -l 
8 7V ' 

(7V 2 -l) (0.0051069297 



(9) 



1 



(7V 2 -l) ( 



128 TV 2 
0.0023152583(50) 



), 



TV 3 

0.002265487(17) 

N + 



0.000794223(19) N). 



Fermionic contributions take the form: 



= 0, 



(iV 2 -l)/l 2 



Nf_ 
N 



(10) 



TV N 2 
4 = {N 2 -l)(h 30 N f + h 3 iJ I + h 32l ±). 

The coefficients h%, /130, /132 depend polyno- 
mially on the clover parameter csw : 



/ l2 =/4 0) +4 1) c S w+4 2) ci' 



(11) 



h,3i=h 



(o)^(i) 



3i 



csw- 



,(2) 



C SW~ 



,(3) .3 



C SW" 



We have calculated the values of h^' , h^' for 



So) u(J) 

l 2 ' "31 

typical values of the bare mass m (see Ref. P for 
detailed numerical values). 



We list below some typical examples of values 
for (E G ), setting N = 3. For N f = we have: 

(l/3).g 2 +0.0339109931(3) 5 4 +0.0137063(2). 9 6 (12) 

For Nf = 2 and m = —0.518106 (corresponding 
to k = (8 + 2m)- 1 = 0.1436): 

c sw = 0.0 : (1/3) g 2 +0.026185200(3) ff 4 
+0.0119649(3) 5 6 , 

c sw = 2.0: (1/3) g 2 +0.013663456(3)# 4 1 ' 
+0.0110200(13)# 6 . 

For N } = 2 and m = 0.038: 

c sw = 0.0 : (1/3) g 2 +0.030438866(3) g 4 
+0.0138181(2) g e , 

c sw = 1-3: (1/3) g 2 +0.025219798(9) g 4 

+0.0129659(5) g e , ( U > 

c sw = 2.0 : (1/3) g 2 +0.01800170(1) g 4 
+0.012948(1) g 6 , 

For N f = 3 and m = 0.038: 

c sw = 0.0 : (1/3) g 2 +0.028702803(5) g 4 
+0.0139032(2) g 6 , 

c sw = 1.3: (1/3) g 2 +0.02087420(1) g 4 

+0.0128495(8) g 6 , ^ Lb > 

c sw = 2.0 : (1/3) g 2 +0.01004706(2) g 4 
+0.013547(1) g 6 . 

3. OVERLAP FERMIONS 

The fermionic action now reads [Hj: 



5 / = EE^N(z,yM 

/ x,y 



(16) 



with: D N = M [l + X(X^X)- 1 ^], and: X = 
Dw — Mo; the sum on / runs over all flavors. 
Here, Z>yv is the massless Wilson-Dirac operator 
with r = 1, and Mq is a free parameter whose 
value must be in the range < Mq < 2, in order 
to guarantee the correct pole structure of Z?n- 

Fermionic vertices are obtained by separating 
the Fourier transform of Dn into a free part (in- 
verse propagator Dq) and an interaction part S 
PHI: 

D N (q,p)=D Q (p)(2Tr) 4 6 4 (q-p) + E(?,p) (17) 

Mq 
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£(g,p)= 



u(p) + w(q) 



[Xi(q,p) 



(18) 



.-^Xo(^t(g,p)Xo( P )] 



w(p) 



- 1 X ( g )Xt(g,p)X (p)] 
w(p)w(g) 



/ 



d 4 k 



1 



(2tt) 4 w(p) + u>(q) w(p) + uj(k) u(q) + w(k) 
-X (q)Xl(q,k)X 1 (k,p) 

-X 1 (q,k)xt(k,p)X (p) 

uj{p) + uj{q) + uj(k) ^ 
w(p)w(g)w(fc) 

X (g)X 1 t (g,fc)X (A ; )X 1 t (fc,p)X (p)l + 0( 5 3 ) 



where Xo,Xi,X2 denote the parts of the Dirac- 
Wilson operator with 0, 1, 2 gluons (of order 
O(g ), Oig 1 ), 0{g 2 ), respectively), and: 

a;(p)=(^si n 2 p M +[^(l-cosp /1 )-M ] 2 ) 1 /2(i9) 



From the form of T,(q,p), we see 
that the 2-gluon vertex splits into 
two different types of contributions: 



Only two fermionic diagrams contribute to 2- 
loops, leading to: 



cf 



0. 



12 N 



(20) 



The quantity hi now depends only on the param- 
eter Mq ; our result is shown in the graph below. 
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It is worth noting that the dependence on Mo 
is smooth all the way to the cndpoint values 
Mo = 0, 2, despite the change in propagator poles 
at these values. 

Typical values of (Ew) '■ Setting N = 3, Nf = 
2, we obtain: 



M = 0.01 
M = 1.99 

N f = 0: 



(1/3) g 1 
(V3) g 2 
(1/3) g 1 



-0.029372693(2) g 4 , 
-0.01975396(10) 5 4 , (21) 
-0.03391099316 <? 4 . 
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